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Abstract — Infiltration is a complex physical process in time 

and space, which is difficult to characterize with precision under 
the intrinsic heterogeneous and dynamic soil conditions. The 
precise determination of the water infiltration in the soil helps to 
minimize the risks of degradation, especially with the processes 
of surface runof and flood. Dashte-Abbas is one of the most 
talented agricultural areas of the Ilam province and Iran 
country, measurement and classification of soil infiltration is 
essential to better manage water resources and the prevention of 
desertification in the Dashte-Abbas plain. Measurement of soil 
infiltration is time and cost consuming and but many models 
have been developed for infiltration estimate, which prevents of 
a lot of time and high costs consuming. Kostiakov is one of the 
models that is being used widely to estimate soil infiltration, 
which as the most widely used model in the planing of soil and 
water. The main objectives of this study were: to; 1) validate the 
Kostiakov model estimates with the measured cumulative 
infiltration and 2) to determine the coefficients of Kostiakov 
model. For this purpose, soil infiltration measured using the 
double-ring infiltrometer in 37 stations with 3 repetitions. Using 
the data obtained, Kostiakov model coefficients were 
determined. Mean values of a and c coefficients were 0.4275 and 
0.7542, respectively. Results showed that the range of variation 
of coefficient of c was large, which demonstrated considerable 
spatial variability in the study area. The normalized root mean 
square error (NRMSE) indicates the total difference between 
the measured and modeled (estimated) values was 
NRMSE=0.15, and performance efficiency was EF= 0.74. Thus, 
it confirms that the Kostiakov model accurately estimates the 
cumulative infiltration in the Dashte-Abbas plain in 
southwestern Iran. 

Index Terms — Valdation model, Soil infiltration, Kostiakov 
model, Double-ring infiltrometer, Dashte-Abbas plain. 

I. Introduction 

Quantifying of soil infiltration as one of the major 
components in the hydrological cycle can be useful in the 
management of catchments. Soils with high restricts, infiltrate 
smaller amount of rainfall into soil and bring about produce of 
the more runoff and flood [1]. Inversly, with the increasing 
infiltrate of water into soil, decrease runoff and flood and to 
reduce the human and financial losses [2]. 

Soil infiltration is a key factor in the rainfall and runoff 

models [3] and an essential factor for increasing agricultural 


production, since an efficient application of water 
fundamentally depends on the infiltration capacity of the soil 
[4]. soil infiltration plays an important role in crop yield and 
leaching of soils in the agricultural aspects [5], [6]. Therefore, 
study and quantifying of soil infiltration is very importance to 
detemine the amount of available water for plant growth, 
additional water needed for leaching and design of irrigation 
systems [7], [8]. 

In order to estimate soil infiltration, many models such as 
Green-Ampt, Kostiakov, Horton and Philip have been 
developed for this purpose. Measurement of soil infiltration is 
time consuming and application of these models prevents of a 
lot of time and high costs consuming [9]. Specifically the 
Kostiakov is one of the models [10] that is being used widely 
to soil infiltration. This model is one of the best models 
because paying attention to all the conditions and factors 
affecting the soil infiltration process [11]. 

Kostiakov model: 

The model of Kostiakov is defined as equation (1) [10], 
where i(t) is the cumulative infiltration (cm ) as a function of 
time, a and c are the equation’s coefficients (a> 0 and 0 <c< 
1). Coefficients of a and c are different and they depond on 
many factors such as soil type, time, ancient moisture, 
hydraulic conductivity of soil and etc [12]. 

m = ct a a) 

Duan et al., (2011) [13] using six different models showed 
Kostiakov and Horton models provided the best predictions 
of cumulative infiltration and average infiltration rate, 
respectively. Zolfaghari et al (2012) [14] eveluation sevene 
infiltration models for estimate cumulative infiltration in the 
four diiferent claases of soil including loam-clayey, 
loam-silty, loam and loam-clay-silty textures through the use 
of ring double. Zolfaghari et al (2012) [14] reported that 
kostiakov and SCS models estimated cumulative infiltration 
properly and imprecise in all the soil classes, respectively. 
Mirzaee et al., (2014) [15] using some models in the different 
areas of Iran, reported that Kostiakov model among other 
models had been the best fitness with measured data. 

Due to low rainfall (220 mm y" 1 ) and its poor distribution in 

the country, Iran is considered among the arid countries which 
faced with water shortage [16]. Around 70% of Iran’s 
agricultural lands are located in arid and semi-arid regions of 
the country. Dashte-Abbas plain with the mean annual rainfall 
235 mm is a arid region of Iran that almost 85% of its 
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agricultural lands have suitable capable for irrigation. Also 
given the important role of soil infiltration in the prevent loss 
of water during irrigation, measuring and quantifying soil 
infiltration to better manage water resources and rising 
irrigation efficiency is essential in this plain. The main 
purpose of this study were: to; 1) validate the Kostiakov 
model estimates with the measured cumulative infiltration and 
2) to determine the coefficients of Kostiakov model as the 
most widely used model in the planing of soil and water, in the 
Dashtabass plain in southwest of Iran. 

II. MATERIALS AND METHODS 
a. Study area 

The Dashte-Abbas plain located about 67km of Dehloran, 
(32’27”N, 47’ 25”E) Ilam province, in southwestern Iran 
(Fig. I). 


Based on data obtained from the meteorological station in 
Dehloran city (southwest Iran), the total annual precipitation 
is 235 mm, the mean annual temperature is 26°C, and the total 
annual open pan evaporation rate in the area are 4300 mm, 
respectively. The climate type is classified as arid region 
according to the De-martonne classification, with a distinct 
dry season during summers and relatively humid during the 
winters. The Dashte-Abbas plain contain Piedmont Alluvial 
Plains soil properties and lies within a relatively flat basin 
physiography. Soils are mainly Entisols based on soil [17] 
with sandy textures (75.7% of sand). Croplands are current 
dominant vegetation type in the plain. The most croplands 
including corn, wheat and barley in this plain. Also almost 
85% of its agricultural lands have suitable capable for 
irrigation (Fig. II). 
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Fig. I. Location map of the Dashte Abbas plain in the southwestern Iran. 


b. Field measurments 

Infiltration was measured using a double-ring infiltrometer 
[18], [19] in the Dashte-Abbas plain. Double-ring 

infiltrometers had been 30 and 60 cm internal and external 


diameter, respectively, and height of 50 cm. Infiltration tests 
were carried out at time intervals of 1, 2, 4, 6, 11, 16, 26, 36, 
56, 76, 106, 136, 166 and 196 min and often have been 
measureds continously untill infiltration rate arrived to a level 
of reletaively constan. determined cumulative infiltration of 
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soil in 37 stations with 3 repetitions and alltogether 111 soil samples were measured. 
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Fig. II. Classification of lands for surface irrigatin in the Dashte-Abbas plain (I class is sutaible for irrigation and VI 


class is not sutaible). 


c. Fitness of Kostiakov model to measured data 

After determined cumulative infiltration in 37 stations (each 
with 3 repetitions) in the Dashte-Abbas plain, we used these 
data to determine coefficients of Kostiakov model. In the first 
defined goal function as equation (2) [20]. Where I(m)j and 
I(p)j are measured data and estimated cumulative infiltration 
by Kostiakov model at j nd time as well as n is the number of 
the paired values, and SSE is sum of square error (cm 2 ). For 
determine coefficients of Kostiakov model, SSE minimum 
amount chosen to the best fit between measured and estimated 
to be created. 


SSE = (2) 

i =1 


d. Validation of Kostiakov model 


2 

measured data including determination factor (R ), 
correlation coefficient (r), normalized root mean square 
error (NRMSE) (Eq. 2) and the performance efficiency 
(EF) were used for model validation (Eq. 3). As described 
the smallest value for NRMSE is zero, indicating that 
there is no difference between measured and estimted 
values. The model’s best performance is at EF=1. Where 
I(m)j and I(p)j are measured data and estimated 
cumulative infiltration by Kostiakov model at jnd time as 
well as I(mean) is mean measured cumulative infiltration 
and n is the number of the paired values. 


NRMSE = 


f X-rrfKm); - I(p)jf 

4 n 


/ (fflSOtt) 



EF = 


EfU - liman) I - If =L (/(p) . -i(mj 

- /(maan)j 



Using the cumulative infiltration data obtained from the 
this study the Kostiakov model was validated. Some 
statistical comparisons between the estimated and 
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III. RESULTS 

a. Determine of coefficients for Kostiakov model in the 

Dashte Abbas plain 

Table I shows the range of variation coefficients of 
Kostiakov model in the Dashte Abbas plain. Mean values of a 
and c coefficients were 0.4275 and 0.7542, respectively (See 
table I) which approximately mean values of coefficients of 
Kostiakov model are high in this plain. Coefficients of c and a 
should be bigger than zero and between zero and one, 
respectively which in this study coefficients were in the range 
mentioned (See table I), precent of coefficient variation 
(CV%) were 31.57 and 51.6 for coefficients of a and c, 
respectively which demonstrated the range of variation of c 
and a coefficients were large in the Dashte-Abbas plain. (See 
table I). Finally, general form of Kostiakov model in the 
Dashte-Abbas plain is as equation (5), which farmers and 
manegers can use it to determine the infiltration in this area. 


i{t) = 0.7641 75 (5) 

b. Validation of Kostiakov model in the Dashta Abbas 

plain 

A significant linear relationship (R = 0.98) was found 
between the measured cumulative infiltration data and the 
estimated lines (Fig. III). Also correlation coefficient (r) was 
0.99. A comparison of modeled (estimated) vs. measured 
cumulative infiltration in the Dashte-Abbas allows an 
overview on Kostiakov’s performance in estimating 
cumulative infiltration (Fig. III). Also as shown (Fig. Ill), the 
measured and modeled (estimated) cumulative infiltration 
values were found near distribution lines. The normalized 
root mean square error (NRMSE) indicates the total 
difference between the measured and modeled (estimated) 
values was NRMSE=0.15, and performance efficiency was 
EF= 0.74 (See table II). 


Table I. Range of coefficients of Kostiakov model in the Dashte-Abbas plain. 



Table II. Quantitative statistical analysis of between measured and modelled (estimated) cumulative infiltration. 




Cumulative infiltration measured {cm) 


Fig. III. Measured and modelled (estimated) cumulative infiltration in the Dashte Abbas plain compare 

with 1:1 line. 

levels of sand (75.7%) in the Dashte-Abbas plain. Also results 

IV . Discussion showed that the range of variation of coefficient of c was large 

^ t t . .... (See table I), which demonstrated considerable spatial 

Result showed that values ot coefficients ot Kostiakov 

t t ,, _ , . , , , _ t , variability in the study area. Measured data were well 

model are high (See table I), which probably due to high 
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correlated with estimated values obtained from the Kostiakov 
model (See table II and Fig. II). Thus, it confirms that the 
Kostiakov model accurately estimates the cumulative 
infiltration in the Dashte-Abbas plain in southwestern Iran. 
These results have been consistent with studies of [13], [14], 
[15]. 
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